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Abstract 

Recent calculations of the order (Za) 4 -jgRy pure recoil correction to hydrogen 
energy levels are critically revised. The origins of errors made in the previous 
works are elucidated. In the framework of a successive approach, we obtain the 
new result for the correction to S levels. It amounts to —16.4 kHz in the ground 
state and —1.9 kHz in the 2S state. 



1 Introduction 



The correction to S levels of hydrogen atom, that is first-order in m/M and fourth- 
order in Za, has become recently a point of controversy. Initially, this correction was 
calculated in Ref. 0. Then, a different result for the same correction was obtained 
in Ref.||. While in both papers it was employed the same (exact in Za) starting 
expression for the pure recoil correction, the methods of calculation, and in particular 
the regularization schemes used were rather different. To resolve the discrepancy between 
two results, in Ref. || an attempt was undertaken to prove the correctness of the earlier 
result of Ref. fl| applying the method of calculation used by the present author in Ref.||. 
An extra contribution due to the peculiarities of the regularization procedure was found 
by the authors of Ref. [||], which exactly compensated the difference of the Ref.@ result 
from that of Ref.|I[. This finding has led the authors of Ref.|| to conclusion "that 
discrepancies between the different results for the correction of order [ZafimjM ) to 
the energy levels of the hydrogenlike ions are resolved and the correction of this order is 
now firmly established" . 

Taking criticism of the Ref. || as completely valid, we nevertheless cannot agree with 
the conclusion cited above. The point is that the authors of Ref.|| emphasizing an 
importance of an explicit regularization of divergent expressions, pay no attention to an 
accurate matching of regularized contributions. 

In fact, one usually starts from an exact expression which can be easily checked 
to have a finite value. Then one has to use different approximations to handle this 
expression at different scales. In this way some auxiliary parameter(s) are introduced 
which enable one to separate applicability domains for different approximations. Finally, 
a necessary condition for the sum of thus calculated contributions to be correct is its 
independence from any scale separating parameter. 

In the present paper we successively pursue this line of reasoning for a recalculation 
of the order (Za) 6 m 2 /M correction to hydrogen energy levels. We discuss only S levels 
since for higher angular momenta levels the result is actually firmly established |4|, 0]. 
As far as the controversy mentioned above concerns details of a regularization at the 
subatomic scale, the result's dependence on a principal quantum number n is also known. 
That's why we perform all the calculations for the ground state and then restore the n 
dependence in the final result. 

To make the presentation self-contained we rederive some known results, using some- 
times new approaches. In Sec. 2 the general outline of the problem is given. Sections 3, 
4 and 5 are devoted to the Coulomb, magnetic, and seagull contributions, respectively. 
The correspondence between various results is discussed in Conclusion. In Appendixes, 
we address a couple of minor computational issues. 

Throughout the paper the Coulomb gauge of electromagnetic potentials and rela- 
tivistic units H = c = 1 are used. Leaving aside the radiative corrections we set Z = 1 
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in what follows. 



2 General Outline 



The first recoil correction to a bound state energy of the relativistic electron in the 
Coulomb field is an average value of the non-local operator, ]5j, || [7], ^ || 

AE rec = --^ J ^((p- D(co,f'))G(r',f\E + u>) (p- D(co,r})) , (1) 

taken over an eigenstate of the Dirac equation in the Coulomb field, 

a 

Hxp(r) = Eip(r), H = ap + (3m- -. (2) 

In (0), p is the electron momentum operator, D(u,r) describes an exchange by the 
transverse (magnetic) quantum, 

-a f d3 k ikf^aak _ _ k(ak) , . 

while 

G(f',f|£; + a;) = \E + u-ap- (3m + ^j 8{r'-r) (4) 

is the Green's function for the Dirac equation in the Coulomb field. The integration 
contour in (0) goes from the minus infinity to zero below the real axis, round zero from 
above and then proceeds to the plus infinity above the real axis. 

As far as we are going to calculate the correction (|1|) perturbatively, i.e. as a power 
series in a, it proves convenient to decompose (|l|) into three parts, 

AE rec = C + M + S, (5) 

namely the Coulomb, magnetic and seagull contributions, corresponding to pp, pD + Dp 
and DD terms from (j^) respectively. 



3 Coulomb Contribution 

It is natural to continuously transform the integration contour into the sum of two 
sub-contours, thus splitting the Coulomb contribution into two terms, 
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where A_ is the projector to the set of negative-energy Dirac- Coulomb eigenstates. The 
former term in (Q) results from the integration along the upper half of the infinite 
circumference and its value is determined by the atomic scale p ~ ma. Being the 
average of the local operator, this term can be easily calculated exactly. The latter term 
in (|6]) arises as an integral along the contour C_, wrapping the half-axis (—00, 0) in the 
counterclockwise direction. To the order we discuss, this term is completely saturated 
by momenta from the relativistic scale p ~ m. That's why it can be calculated without 
any regularization [1], [J: 

I w? eft 

~ m = —tf-- (7) 

4 Magnetic Contribution 

Using the identity 

(pGD + DGp) = 1([P, H]GD + DG[H, p\ + {p, 3}) , (8) 

which follows directly from the equation for the Green's function, we can extract from 
the general expression for the magnetic contribution, 

its local part, 




Due to the rapid convergency of the integral in (|9D at the infinity, the integration contour 
can be reduced to C_. By virtue of the virial relations (see || and references therein), 
the sum of local parts of the Coulomb and magnetic contributions takes a simple form 

Physically, this contribution to the recoil correction is induced by an instantaneous part 
of the electron-nucleus interaction. 

4.1 Long Distances 

Immediate integration with respect to u in (§) gives ||: 

M = -— f ^ /n/V |m)H V |m)(m| \ A ™ k r i &\ (12) 
M J (2ir) 3 \ \+ k + E m — E ^ E - E m + k ) k /' K ' 

where S+r-) stands for the sum over discrete levels supplied by the integral over positive- 
negative-) energy part of the continuous spectrum. 
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4.1.1 Positive Energies 



In the leading nonrelativistic approximation, the first term in Eq.([E|) reads, 

^^/(S( ?e(f ' ri£ - fc, ^>' (i3) 

where Q [r ', r\E — k) is the Green's function for the Schrodinger equation in the Coulomb 
field, and the average is taken now over the nonrelativistic wavefunction. For the ground 
state, we work with 



^(f) = ^^e— , E = -^f. (14) 

Only p-wave term from the partial expansion, 

£(f',r» = y £{-) l (2l + l)Pi(m)g l (r',r\u), (15) 
i 

survives the integration over the angles: 

M+ = --^ dkk dx(l-x 2 )(£i(r',r|£-A;)e fc y (16) 
Mn Jo J-i ^ 1 

For the nonrelativistic Green's function in the Coulomb field we use the integral repre- 



sentation from the paper [TO 
Qi r'.r 



--) 


im r 


2m J 


2it\fr J r Jo 



: J 3\ — : • U'J 

sin s 1 — e » V sin s j 



The integrals over r and r' in ( |Tq) are easily calculated after expanding the Bessel 
function into the power series. The result can be expressed in the form, 

2 7 3m 5 a 6 r°° dk k r 1 , . 9x r , t l ~ ma / K 1 , , 

M+ = — — / — =- / dx(l-x) / dt- — tt-t^t, (18) 

+ Mvr Jo k J-i v ' Jc {a-btyi- e 2m ~ V ; 



where k = ^2m{k — E), the contour C is the unit circumference \t\ = 1 directed clock- 
wise, and 

/ ma\ ( ma ikx\ , / ma\ ( ma ikx\ 
a= 1 + 1 + , b= 1 1 + — • 

V K J \ K K I \ K / \ K K I 



Integration by parts conveniently extracts from the last integral in ([18|) the terms non- 
vanishing at large momenta: 

M + = — - / dy(l-y 2 ) dx(l-x 2 )F(x,y), (19) 



Mvr Jo 

where 



{,y} b{a-b) 3 b 2 {a-b) 2 ab\a - b) aW a% Jo 1 - H' K ' 
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and the new integration variable y = ma/ k is introduced. Since 



*=fl^, (21) 
k 2 y 



to get a power series expansion of ( ]T9| ) with respect to a up to the first order, we need 
an expansion of the integrand with respect to y also up to the first order (note that 
a — b = Ay — 2ikx/ ' k): 

(1 - »")*■(,, y) „ ^ - ^ + i - £ + yMa _ b) . (22) 

Here the last term emerges as a result of expansion of the integral in (|20D , 

1 dtt- y 1 / fc\ 

^-S7= i F l,l-y;2-y;- , (23) 

o 1- ^ 1 -y \ a J 

where F(l, 1 — y; 2 — y; b/a) is the Gauss hypergeometric function. Integrating now ( |22|) 
first with respect to x, and then with respect to y from to some y (a 1 / 2 <^ y <^ 1), 
we obtain, 

f° rfo(l - y 2 ) /' - x 2 )F(x, y )«^--L-ll n M + _L (24) 
Jo J-i 32a 48t/o 12 a 48 

1 2y 2yl , , 3un tto; 

— + — + — ln4y - — ■ 

9 3 3 y 4 32 

On the other hand, we can neglect a in F(x, y) on the interval [yo, 1]. In the sum of two 
integrals, the dependence on the auxiliary parameter y disappears, and we come to the 
result, 

m 2 a 5 \ 7T 8 , 1 8 , Ry 16 , 32 ] 



Here the Bethe logarithm is introduced |TT 



'I dyy a + vm-v) = ln (EA7, 



16/ dyy^ ^ ' = i n _^ + 2 ln2 + ^. (26) 

lo (1 + y) 4 (l - y) 6 

In (|25|), the order a 4 term is just the lowest-order contribution to (PH), the order a 5 
terms are in accord with the result of Salpeter [fT^| , while the order a 6 term coincides 
with the retardation correction, found in ||,Eq.(14) by the different method. 

It can be easily seen that the order a 6 contribution to the positive-energy part of 
(|12|) is exhausted by the sum of those to ( |i~0| ) and (^). Actually, relativistic corrections 
are at least of the a 2 relative order. The effect of retardation reveals itself starting from 
the a 5 order (^). Hence the relativistic corrections to the retardation are at least of 
the a 7 order. 



5 



4.1.2 Negative Energies 



Virtual transitions into negative-energy states give rise to the second term in (|12|). In 
the leading nonrelativistic approximation, it equals [0, 

Am 2 M J (2vr) 3 \k' 2 k 2 /' V 1 

where = k — k'{kk')/k' 2 , k' = p' — p — k, p and p' being the arguments of the 
wavefunction and its conjugate respectively. The integral over k diverges logarithmi- 
cally (leading linear divergency vanishes due to the numerator which at k—>oo becomes 
transverse to itself, and hence rises only like k, not k 2 ). To treat this divergency we 
use the following formal trick 0: subtract from ( p7[) the same expression with k' 2 + A 2 
substituted in place of k' 2 . For A ^> ma, the subtracted term is completely determined 
by a scale much less than the atomic one, so that we will find that term below using a 
relativistic approach. 



The regularized version of fl2T|) can be written in the form 

a 2 I . , . r d 3 k Arak j ( . k'M\ /4tt 4tt , > 

In the coordinate representation, the integral above is 



After substitution into (p8|) it gives 

Finally, the result of trivial calculation of the average over the ground state reads, 

m a /„, e 



where e = X/2m. 



M--M r . = -^r- 2 In - - 1 , (31) 
M \ a 



4.2 Short Distances 

Since in the nonrelativistic approximation the subtracted term, A4 r _, is ultraviolet di- 
vergent, we have to calculate it beyond this approximation, i. e. using a relativistic 
approach. It proves more convenient in this approach to postpone the integration over 
uj to the last stage of calculation. As we will see below, the reversed order of integration 
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(first over space variables, then over frequency) makes the calculations quite simple. The 
fee for the technical advantage is that a regulator contribution is calculated not only for 
the negative-, but for the positive-energy part of M also. Surely, the instantaneous 
contribution can be left aside, so that only two first terms from the r.h.s. of (§) are 
considered below. 

For the subtracted term, we have the new expansion parameter, ma/A, and hence the 
Coulomb interaction during the single magnetic exchange can be treated perturbatively. 
The order ma 6 /M contributions arise due to only two first terms of the Green's function 
expansion in the Coulomb interaction, and . Let us begin with the second 
contribution: 

M'a = li £i(Kfl]G«a r ). (32) 



M Jc- 2ni u 

Here 



*r 



D = / — — -e 



d 3 k 47ra5i 



(2vr) 3 k 2 + \ 2 -uj 2 ' 
and we can neglect atomic momenta in comparison with A and m: 

a 3 ip 2 f dcu I Anp ' 2m + uj + ap' An u + ap 4ira p \ , . 

= tvM~ Jc- ~ioo P' 2 - fi 2 g 2 " p 2 -W p 2 -IC 2 /' ^ ' 

The notations of || are used: ip 2 = |^(0)| 2 , the angle brackets denote here integrations 
over p and p' together with the average over the spinor u a — 8 a \\ q — p' — p; and 



K = Vio 2 - A 2 , Q = V2mu + J 2 . 
The average over the spin degrees of freedom gives 

((2m + u + ap'){uj + ap)a p p') = up' 2 p = ojp' p q. (34) 
Then, after transition to the coordinate representation we get 



2a V f du f*> j f Q e lUr - 1 



Mr = / — / dr\ di n,- 

G mM Jc- iu Jo \ tt 2 r 1 3 



n 2 

The integration over r is simple but lengthy. It results in 



d d \ e iQr - 1 



(35) 



Ma = -P^l -|^-Slnfl + ^ + (^^) + 21nfl + ^). (36) 
G 2mMJc-iu\VL Q 2 \ KJ v ; V ttj f K J 

Here the contour of integration goes counterclockwise around the cut connecting points 
—2m and —A. According to the Feynman rules, Q = while JC = +(— )\JC\ on the 
lower (upper) edge of this cut. Since the integrand is regular at small ui, we can put 
A = (recall that A<Cm) and get 

a 3 ^ 2 f 1 ( y/l — x 1 / x 1 \ 3 na 3 ip 2 



Mr = — tt I dx — -arctan./ . = — — . (37) 

' ' mM Jo \ x 3 / 2 x 2 V 1 - x l x (i _ x ) } 2 mM 
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To calculate a contribution due to we have to account properly for the wave- 
function's short- distance behavior: 

. , r a 3 ip 2 ( du I I Anp' 2m + uj + ap' Ana q 

* = !c. ^ \\]/* p' 2 - Q 2 q 2 -K 2 [ ^'' 

An dpi uj + dp 1 Anq\ 2m + dp Arc \ 



p ,2 _ ^2 p ,2 _Q2 q 2 J p 2 p 2 

Averaging over the spin part of the wavefunction, we obtain 



_ «V (du I (Am \ An An Anp\ 



An An Anp 2 , 
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( P ' 2 -ic 2 )(p' 2 -n 2 ) q 2 p 

Again, the six- dimensional integral over p and p' turns into a simple integral over r in 
the coordinate representation, and equals 



2a 3 ip 2 [du [(Am ATI, / Sl\ 1 / /C 
= I — \ ( + 1) — In ( 1 + - ) + — In [ 1 + - 



M J i I V u 



tt 2 V K) K, 2 \ ttj fi/C 

1 ln5 



(40) 



Finally, the integration along the same contour as above gives for non-vanishing in the 
limit e— >0 terms: 

kjr m 2 a 6 (2 32 f°° d6 n . l\ 



We see that as expected the logarithmic in e term cancels the corresponding one in (pl|). 
The more singular in e terms can only be the result of the regularization procedure 
applied to the positive-energy contribution (|25|). As far as the latter is non-singular at 
short distances, this procedure is actually unnecessary, i. e. it can produce only positive 
powers of ma/ A. An explicit calculation can be found in Appendix A. 



4.3 Total Magnetic Contribution 

So, in the sum of all contributions due to a single magnetic exchange, any dependence 
on the scale separating parameter e cancels away, and we get 

— (-i+^-i- §)■ <-> 

Here —1 in the r.h.s. is due to the (long-distance) effect of retardation (see ( p5|) and 0, 
Eq.(14)), 2 In ^ — 1 comes from the whole range of scales from ma to m, while —3/2 is 
the short-distance contribution. 
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5 Seagull Contribution 



5.1 Long Distances 

Again, the best suited way to analyze the atomic scale contribution is to begin from 
taking the integral with respect to to. It proves that in the order of interest, only 
positive-energy intermediate states are to be considered ||: 

a 2 r d 3 k 1 Att 2p' k , + ia x k' Arc 2p k + ia x k\ 
+ = 2Mi (2^\F 2^1 F 2^ /' { ' 

A simple power counting shows that only bilinear in k and k' term gives rise to the 
ultraviolet divergency. To regularize this divergency, we subtract from the divergent 
term the regulator contribution, which at large distances equals to 

a 2 I Auk' Auk \ , . 

(44) 



Am?M \ k' 2 + A' 2 k 2 + X 2 / ' 

while ma <C A, A' <C m. In the coordinate representation, the regularized version of (f43"D 
is 

5 + " 5 ' + = 4^M ( 2 4 ?+ h ~ ( V£ ^) ( V V 1 ) ) ' (45) 
The average over the ground state reads (e' = X' /2m): 

m 2 a e \ e' 2 + e'e + e 2 , e' + e 2e'e ) . , 

S+-S r ,=— — (2 — + l-21n + - -}. 46 

+ + M \ a(e' + e) a (e' + e) 2 J v ; 

Here 1 appears due to the non-singular operator pr~ 2 p. The first term in the curly 
brackets represents the regulator contribution to the previous order. In Appendix B, an 
appearance of this term as a short-distance contribution to the ma 5 /M order is shown 
explicitly. In what follows we calculate the subtracted term, whose non-relativistic 
version (H) is ultraviolet divergent, in the framework of a relativistic approach. 



5.2 Short Distances 

Just like in case of the single magnetic exchange, only two first terms of the Green's 
function expansion in the Coulomb interaction contribute to the m 2 a e /M order. For 
the GW's contribution we have, 

r a 3 ip 2 f duj I Ana p i u + ap' An u + ap Ana p \ 

G = 2vrM Jc. ~T\p' 2 - K' 2 p' 2 - Q 2 If p 2 - Q? p 2 - K 2 ) ' ( ' 

Calculation along the same lines as in the case of Ai r G gives the result 

S «= M^ (41ll2 - 2 >' (48 > 
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which is non-singular in the limit A, A'— >0. 

As for the contribution due to G^°\ it can be extracted from 

a 3 ip 2 r du I 4:7idp> uj + dp' Amd q 2m + dp \ , 

2nM Jc. T \p>2-lC' 2 p' 2 - Q 2 q 2 -K? [f + 7 2 ) 2 / + ( ~ } ' ( } 

as a zeroth-order term of the Laurent series in 7 = ma (this series begins with an 
order I/7 term describing the seagull contribution to m 2 a 5 /M order at short distances 
discussed in Appendix B). The average over the spin part of the wavefunction is 

(2mud p ,d q + d p ,dp'd q dp) = -(u 2 + \p a - n 2 }) (1 + + W (50) 



The term in the square brackets can be omitted. In fact, the corresponding part of (f49|) 
does not depend on m and hence (merely on dimensional grounds) contribute to the 
m 2 a 5 /M order only. Then, the first term gives the non-singular in the limit A, A'— >0 
contribution: 

^^ (1 _ 41n2) . (51) 

\ p' 2 q 2 J Mm K J ( J 

Finally, analysis of the last term from (|50|) deserves more care since here we have the 
infrared singularity. Being integrated over the space variables, this term gives: 

2p'q^ 2 -^( — {m,K)-f{K\K)}, (52) 
Mm JC- iuj 

where 

/(*,»)= in (l + I) (53) 

(recall that K,' = \J uj 2 — A' 2 ). For e < 1 we obtain 

2a 3 ib 2 f du „,„ ira 3 ib 2 ( , 1 \ . . 

' ./(n j /C) = — — -21n- + 41n2-l . (54) 



Mm JC- iuj Mm \ e 

Calculation of the integral with /(/C',/C) is a bit more cumbersome since it does not 
contain a small parameter. The contour C_ for this integral encompasses in the counter- 
clockwise direction the cut connecting the points —A and —A'. Continuous deformation 
of C_ leads to the following equation: 



du ... = / du . . . — 2ni Res ... — 2ni Res . . . , (55) 

C- Jc+ ^=0 ^=00 

where ... stands for f(K.',JC)/u, and the contour C + goes in the clockwise direction 
around the cut connecting the points A and A'. Using the evident relations, 

du . . . = — du . . . , (56) 

C+ JC- 

Res !/(£', K) = /(A', A), (57) 

U=0 uj 

Res -/(£', K) = -ln2+-, (58) 
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we come to the result 

2o?^ 2 f duo . J . ixo?i) 2 ( 2s 2ee' V 



5.3 Total Seagull Contribution 

As can be seen from (|46"|), fl48|), (|51~D, (|54"D, and (|59|), the total seagull contribution to 
the m 2 a 6 /M order does not depend on the scale separating parameters A and A', and 
equals 

S a , = lllfl (l-21n- + - + 41n2-2) , (60) 



M V a 2 

where 1 comes from the long distances, 4 In 2 — 2 from the short ones, while the remaining 
terms gain their values on the whole range of scales, from ma to m. 



6 Conclusions 



In complete accord with the result of Ref.|13[, the total correction of the rrror/M order 



does not contain In a. It consists of two terms 

m 2 — E 



2M 



ma 



+ w' 2 " l2 - 3 >- (61 > 



The former term is completely determined by the atomic scale and depends non-trivially 
on a principal quantum number n, 

m 2 - E 2 m 2 a % (\ 3 2 1 \ 

(62) 



2M 



2Mn 3 V4 An n 2 n 3 



As for the latter one, our calculations show that it has its origin at the scale of the order 
of m. 

The correction ([H]) shifts the hydrogen ground state by —16.4 kHz and 2S state by 
— 1.9 kHz. These figures are well comparable with the uncertainties of the recent Lamb 



shift measurements 114 



The result (|6~i"l) differs from those obtained in Ref.fll], |j and in Ref.@. Let us first 
discuss the origin of the difference in the latter case. In 0, it was erroneously assumed 
that the cancellation of singular operators at the atomic scale does not leave a non- 
vanishing remainder. The present calculation shows that due to a difference in details 
of a cut-off procedure used to regularize the average values of singular operators, some 
finite contributions do survive the cancellation process. 

Unfortunately, the same error was repeated in Ref. |§. The long-distance contribution 
was found there in the framework of some particular regularization scheme. Then it was 
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added to the short- distance contribution calculated in Refs. [|I], Q] by completely different 
regularization procedures. The regularization dependence of the results obtained in || 
can be seen, for example, in Eq.(29) of Ref.||, where the integration over k' being limited 
above by a parameter a' gives rise to a finite (depending on a' jo) contribution to the 
result. 

The error made in Ref . Q is a computational one. It is caused by inaccurate treatment 
of the frequency dependence in the integral (42) of Ref.|]l| (ironically, by an evident 
typographical error, just the important factors (uj 2 — k\)~ l and (uj 2 — k 2 )' 1 are skipped in 
Eq.(42) of Ref.0). In what follows we rederive the result of the present work employing 
the regularization scheme used by the authors of Ref. 

First of all, the result for the long-distance contribution (46) of Ref. |I[ ("the third 
term") is in accord with the result of the present work ( 1 in ([46])). 

As for the remaining contributions, let us begin with one general note. In their 
analysis of the integral (42), the authors of Ref.|l| use the symmetrization in oj, since, 
as they wrote, "generally there are three regions of photon energy uj ~ a 2 , uj ~ a, 
and oj ~ 1 that give a contribution and this middle region is almost eliminated by the 
symmetrization". In order not to discuss here whether the middle region is eliminated 
or not, we would like to recalculate the contributions of the first and the second terms in 
(43) of Ref. [Q] without the symmetrization in oj. As far as the symmetrization procedure 
is no more than a technical trick, a result of calculation should not depend on whether 
this procedure is applied or not. 

To get the high energy part of the first and second term contribution, we put e' = 
e = in ([S|) and cut off the low energy end |u;| < me from the contour C_. Then the 
result for the short- distance (high energy) contribution to the integral (42) of Ref. [0] can 
be obtained: 

2 fi 

ma e , , 

AE = — — 2 In-. 63 
M 2 v ; 

The sum of the order m 2 a 6 /M contributions to Eqs.(51), (54) and (57) of Ref.|l| is two 
times smaller. An extra factor one half emerges there due to the symmetrization in uj, 
since the contribution of the contour C+, wrapping the half-axis (me, oo), vanishes. 

Turn now to the low energies. Only the second term of Eq.(43), Ref.0, contributes 
there. According to (42) and (43) of Ref. [0], this contribution (with typos corrected) is, 

A£ = ®— f r^h. ffk (1L m 

Mm Jc L 2m J (2tt) 3 J (2tt) 3 J (2vr) 3 1 } 

t 471"^! 1 47rA; 2 r> 

WVP + kl)j2 2 O 2 72 i^U 3 + k V- 

k{ — uj 1 2muj — p 2 ^2 — uj 2 

Here the contour Cl goes from —me to below and then from to me above the real 
axis. Recall now that the high energy contribution (|5"3"| ) is calculated on the assumption 
that e>a. It means that in (|64]) we can neglect p 2 which is of the order of (ma) 2 , in 
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comparison with 2muo which will be shown below to be of the order of m 2 a. Then we 
can easily come to the coordinate representation and get, 

AE- ^ I dU 1 /fv e<Hr V *\ (65) 
~ 2Mm 2 Jc L 2wl uo - \ \ ~V~ ) ~ / ' ^ ' 

Since the integration contour does not wrap the zero point, we can safely add the operator 
— 1/r 4 which is annihilated by the uo integration. The result of taking the average over 
the ground state is 

2m 2 a 6 r duo 1 ( ( i\uo\\ 2i\u\ 1 3 / uo \ 2 1 

AE = — — / 2 In 1-— tt + -\ n — r 

M Jc L 2m u - \ \ ma J ma 1 - — 2\maJ 1 - *H 



ma 



(66) 

Here we see that the natural scale for uo is in fact ma. Since |u;| is positive on the lower 
half of Cl, the integral above written in dimensionless units reads, 



m 2 a 6 f e / a , A 4 x 2 » 

Ai? = — — / ax — arctanx — 3 . (67) 

ttM Jo \x 1 + x 2 1 + x 2 I 



The result of integration, 



m 2 a 6 ( e , e 1\ 
AS=— — (3 21n- + - , 68 



being added to all the other seagull contributions, gives for the order m 2 a 6 /M seagull 
correction: 

m 2 a 6 / , 1 1\ . , 

5 „« = _(_ 21n _ + 21ll2 __), (69) 



in complete accord with the result (60) of the present work. 
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Appendix A 

Extra terms in ( f4T|) should be canceled by the regulator counterpart of ( |T3"D which 
differs from ([13|) by \Jk 2 + A 2 placed instead of k. Just like in the main text we approx- 
imate the sum over positive-energy intermediate states by the nonrelativistic Green's 
function and the matrix element of a by p/m. Within this approximation, the regulator 
contribution is 

M + = -Llwr (*"■ ^ - ^) vSSift) • (70) 



13 



After the transformations, the regulator version of the expression (|T^) is 



2 7 3m 5 a 6 r° dk k 2 f 1 , , Qs r , t 1 "™ /* 

.Ml = 



Mtt 



Jo ^L dx{1 - X) Jc d \a-btri-e^ > (71) 

where k = ^2m{u) — E), to = y/k 2 + A 2 , the contour C and functions a and b are defined 
in the text. Only singular terms of the expansion (|22| ) operate at distances of the order 
of A -1 . For those terms the integrals over k and x become elementary and give 

kjr m 2 a 6 ( 2 / e \ 2 32 r°° d6 } 



So, the second and the third terms coincide with the corresponding terms in (|41|). The 
new singularity ~ e~ 2 is nothing but the regulator contribution to the instantaneous 
part of the magnetic exchange ([R]) : 

1 /r^ftr X \ /p'i+Pi Pj+Pj S ij ~q i q J /q 2 \ 

— ({p,D (0,r)}) ~ 4^—-^- q2 + x2 J (73) 

ma /, £ 



= 2 ln-- 

Me 2 V a 

Appendix B 

The leading contribution to ( |49|) is 

8m 3 a 5 r du r d 3 p 47r 47r 
M ic_ 2ri J (2tt) 3 p 2 - Kl 2 p 2 -ft 2 p 2 -K?' ' 

After the integration with respect to p it turns into 

or m « 5 1 /" , ( 1 1 \ ,nr\ 

s = m^^Jc. duJ u [nTlc ~ nTlc) ■ (75) 

Up to terms of the first order in e, e', one gets 

m 2 ^ / 2 g 2 ln(2/ g )-^ 2 ln(2/ £ Q _ ^ + + £ 2 \ 
Mtt V e 2 -e' 2 e' + e / ' 1 J 

The last term compensates the leading contribution to (|46|). 
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